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ABSTRACT 

By combining the data of the Two Micron All Sky Survey, the Wide Field Infrared Survey 
Explorer and the AKARI satellite, we study the infrared colour properties of a sample of 2712 
nearby radio-luminous galaxies (RLGs). These RLGs are divided into radio-loud (RL) active 
galactic nuclei (AGNs), mainly occurring at redshifts of 0.05 < z < 0.3 and star-forming- 
dominated RLGs (SFGs), mainly occurring at redshifts of 0.01 < z < 0.15. RL AGNs 
and SFGs are separately distributed in the ([3.4] —[4.6])—([4.6]—[12]) two-colour diagram, 
in which the RL AGNs display a double-core distribution, and the SFGs display a single¬ 
core distribution. SFGs have a redder [4.6] —[12] colour than RL AGNs due to the significant 
contribution from the dust component of SFGs. We find simple criteria of MIR colour sep¬ 
aration between RL AGNs and SFGs such that: 95% of RL AGNs have [4.6] —[12] < 3.0 
and 94% of SFGs have [4.6] —[12] > 3.0. We also analyse the MIR colours of RL AGNs 
divided into low- and high-excitation radio galaxies (LERGs and HERGs, respectively). 
The ([3.4] —[4.6])—([4.6]—[12]) diagram clearly shows separate distributions of LERGs and 
HERGs and a region of overlap, which suggests that LERGs and HERGs have different 
MIR properties. LERGs are responsible for the double-core distribution of RL AGNs on the 
([3.4] —[4.6])—([4.6] —[12]) diagram. In addition, we also suggest 90—140)xm band spectral 
index Qf(90,140) < —1.4 as a criterion of selecting nearby active galaxies with non-thermal 
emissions at far-infrared wavelengths. 
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1 INTRODUCTION 

Radio-luminous galaxies (RLGs) form a subclass of active galax¬ 
ies. The radio emissions from galaxies are thought to originate from 
star formation and the accretion process around supermassive black 
holes in active galactic nuclei (AGNs). AGNs emit photons over 
the entire electromagnetic spectrum and the radio emissions from 
AGNs are produced by synchrotron radiation of the relativistic jets 
and the lobes formed by interaction of the jets with the surround¬ 
ing media (e.g. Antonucci 1984; Urry & Padovani 1995; Cleary et 
al. 2007). The radio emissions from the star-forming region are be¬ 
lieved to originate from free—free emission in the ionized hydrogen 
(H II) region and the synchrotron radiation by relativistic electrons 
in supernova remnants (Harwit & Pacini 1975). Most active galax¬ 
ies may contain two components: AGNs and star-forming regions. 
For example, NGC 1068 shows an AGN at the centre surrounded 
by a ring-shaped star-forming region (Telesco et al. 1984). 

Both AGNs and star-forming regions require a large supply of 
dusts and gases as fuel. AGNs have different observed subclasses 
that can be explained by a unification model containing a super- 
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massive black hole surrounded by an accretion disc with jets and a 
torus of dust and gas (Antonucci 1993). The dusty torus not only 
supplies the fuel of the accretion disc, but also reprocesses X-ray 
emissions from a hot corona above the accretion disc and then emits 
infrared radiation. Therefore, the infrared bands are important for 
understanding the accretion environment of AGNs. In star-forming 
regions, dust and gas clouds are heated by young massive stars and 
also emit mid-infrared (MIR) and far-infrared (FIR) radiation. 

Active galaxies display various properties at various infrared 
wavelengths. At near-infrared (NIR) wavelengths, emission of 
AGNs is significant, although the contribution of host galaxies to 
the spectral energy distribution (SED) cannot be neglected. For ex¬ 
ample, the NIR spectra of quasars are roughly represented by a 
power-law (PL) continuum superposed by a composite blackbody 
(BB) spectrum of the stellar population (Elvis et al. 1994; Richards 
et al. 2006; Park et al. 2010). The composite BB spectrum peaks at 
~ 1.6/rm. Therefore, NIR colour is also used to distinguish AGNs 
from normal stars (Kouzuma & Yamaoka 2010). At MIR wave¬ 
lengths, dusts in H II regions (or tori), polycyclic aromatic hydro¬ 
carbons (PAHs) in photodissociation regions and AGNs comprise 
potential sources of MIR emission of active galaxies (e.g. Lacy et 
al. 2004; Ichikawa et al. 2014). At EIR wavelengths, it is reported to 
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observe cooler 30 K) interstellar regions heated by stars in host 
galaxies or the outer part of the torus in the AGN (e.g. Ichikawa et 
al. 2012). 

In the present study, we summarize previously reported NIR, 
MIR and FIR properties of radio galaxies. Strong infrared emission 
lines have been commonly detected in high-redshift radio galax¬ 
ies (Rawlings et al. 1989; Bales & Rawlings 1993; Iwamuro et 
al. 1996; Evans 1998). For example, a low-resolution NIR spec¬ 
trum (1.14—1.70/rm) of 4C 40.36 {z = 2.269) shows distinctively 
strong [O III]A4959/5007 doublet emission lines that significantly 
contrihute to the Ff-band image (Iwamuro et al. 1996). A concen¬ 
tric distribution of strong oxygen line-emitting clouds is suggested 
by the 7F-band image of 4C 40.36 (Iwamuro et al. 1996). 

Benefiting from observations implemented by the Infrared 
Space Obserx’atory (ISO) and the Spitzer Space Telescope , many 
previous reports have focused on the MIR properties of high- 
redshift powerful radio galaxies. Siehenmorgen et al. (2004) used 
ISOCAM data recorded on hoard the ISO to construct the SED of 
3CR radio galaxies at redshifts of z^2.5. They determined that the 
MIR properties of 3CR radio sources may include synchrotron ra¬ 
diation from AGN, stars of the host galaxy or dust. Cleary et al. 
(2007) measured the MIR radiation from a sample of extremely 
powerful 3CRR radio galaxies at redshifts of 0.4^zS(1.2 based 
on Spitzer observations. They hypothesized that Eanaroff—Riley 
type II quasars are more luminous than galaxies in the MIR range 
because of a combination of Doppler-boosted synchrotron emis¬ 
sions in quasars and dust extinction in galaxies, both of which are 
orientation-dependent effects. De Breuck et al. (2010) presented a 
Spitzer survey of 70 radio galaxies at redshifts of l<z<5.2 and 
measured the rest-frame 1.6 ^m emission from the stellar popula¬ 
tion and the hot dust emission associated with the active nucleus. 
Leipski et al. (2010) obtained rest-frame 9—16 pm spectra of 11 
quasars and 9 radio galaxies from the 3CRR catalogue at redshifts 
of l<z<1.4 based on Spitzer observations. They determined that 
the mean radio-galaxy spectrum shows a silicate absorption fea¬ 
ture, whereas the mean quasar spectrum shows a silicate emission 
feature (Leipski et al. 2010). 

Based on Infrared Astronomical Satellite (IRAS) data, 3CR 
radio sources were measured at MIR and EIR wavelengths. The 
results confirm that MIR and EIR emissions from 3CR quasars 
and radio galaxies are anisotropic and that the non-thermal beamed 
emission of the core component partly contributes to the FIR emis¬ 
sion of the 3CR radio sources (Heckman et al. 1992; Hes et al. 
1995). Dicken et al. (2009) presented the deep MIR (24 pm) to 
FIR (70 pm) Spitzer photometric observations of a 2Jy sample of 
powerful radio galaxies (0.05<z<0.7). They determined that di¬ 
rect AGN heating is energetically feasible for the dust producing 
both MIR and FIR continua (Dicken et al. 2009). Hardcastle et al. 
(2010) used //erac/reZ-Astrophysical Terahertz Large Area Survey 
data to analyse the FIR properties of radio galaxies. They reported 
that these properties are similar to those of a comparison popula¬ 
tion of radio-quiet galaxies matched in redshift and Ai-band abso¬ 
lute magnitude. 

However, little attention has been paid to the infrared colour of 
nearby RLGs, the activities of which are dominated by either AGNs 
or star formation (Best et al. 2005; Best & Heckman 2012). This 
paper will focus on the infrared colour properties of nearby RLGs, 
including that reported by Best et al. (2005) based on sample of 
2712 RLGs. The main objective of this study is to examine the dif¬ 
ferences in infrared properties of AGN-dominated RLGs and star¬ 
forming region-dominated RLGs by using this large known sample. 
Additionally, AGN-dominated RLGs include two different activ- 
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Figure 1. Redshift distribution of all sources studied. Black line represents 
radio-loud AGNs (RL AGNs); red line represents star-forming galaxies 
(SFGs). 

(A color version of this figure is available in the online journal.) 


ity modes (Hardcastle et al. 2007; Best &. Heckman 2012): high- 
excitation radio galaxies (HERGs) and low-excitation radio galax¬ 
ies (LERGs). HERGs are powered by radiatively efficient standard 
thin discs at high accretion rates, whereas LERGs are fuelled at 
low accretion rates by radiatively inefficient accretion flows, i.e. 
advection-dominated accretion flows (ADAEs; Narayan & Yi 1995; 
Hardcastle et al. 2007; Best &. Heckman 2012; Yuan & Narayan 
2014). Giirkan et al. (2014) analysed the MIR properties of HERGs 
and LERGs by using four complete samples (3CRR, 2Jy, 6CE and 
7CE). They pointed out that HERGs are separated from LERGs in 
the MIR-radio plane (Giirkan et al. 2014). We further analyse the 
MIR colours of HERGs and LERGs by using a sample of nearby 
RLGs. 

Throughout this paper, we adopt a flat universe with a Hubble 
constant Ho = 70 km s“^Mpc~^, Ha = 0.7 and Hm = 0.3. 

The paper is organized as follows: Section 2 describes our 
data; Sections 3—5 discuss the colour properties; and Section 6 pro¬ 
vides a summary. 


2 DATA 
2.1 Sample 

All sources studied here were selected from a catalogue of 2712 
RLGs (Best et al. 2005). This large sample of RLGs was derived 
by Best et al. (2005), who used the second data release of the Sloan 
Digital Sky Survey (York et al. 2000; Stoughton et al. 2002) to 
cross-correlate with combined data of the National Radio Astron¬ 
omy Observatories Very Large Array Sky Survey (Condon et al. 
1998) and the Eaint Images of the Radio Sky at 20 cm (Becker et 
al. 1995) survey. Both radio surveys achieved higher angular reso¬ 
lutions than those recorded by their predecessors and covered sub¬ 
stantial fractions of the sky, down to flux densities at the millijansky 
scale. To remove the very extended galaxies. Best et al. (2005) ex¬ 
cluded galaxies with redshifts of z < 0.01; the 2712 RLGs are 
located at redshifts of 0.01 < 2 < 0.56. These samples are from 
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Table 1. Infrared photometric data of nearby radio-luminous galaxies (RLGs) 


RA(2000) Dec.(2000) 
(1) 

(2) 

E(B - V) 
(3) 

J (mag) 

(4) 

H (mag) 

(5) 

Ks (mag) 

(6) 

W1 (mag) 

(7) 

W 2 (mag) 

(8) 

W3 (mag) 

(9) 

W4 (mag) 
(10) 

90pm ay) 

(11) 

14:0pm ay) 

(12) 

146.956! -0.3423 

0.1347 

0.1343 

15.882(0.109) 

15.120(0.160) 

14.526(0.119) 

14.054(0.035) 

13.829(0.044) 

12.040(0.000) 

8.452(0.000) 


/ 

146.1436 -0.7416 

0.2039 

0.0606 

15.654(0.115) 

14.982(0.133) 

14.416(0.135) 

13.477(0.025) 

13.217(0.030) 

12.156(0.000) 

9.056(0.000) 


/ 

146.7371 -0.2522 

0.1305 

0.1327 

15.482(0.094) 

14.798(0.104) 

14.209(0.093) 

13.505(0.026) 

13.320(0.032) 

12.031(0.274) 

8.996(0.000) 


/ 

146.3738 -0.3684 

0.0529 

0.1118 

15.137(0.094) 

14.474(0.130) 

13.948(0.087) 

12.655(0.024) 

12.492(0.025) 

8.622(0.029) 

6.205(0.068) 


/ 

145.6012-0.0014 

0.1459 

0.0625 

15.958(0.115) 

14.970(0.099) 

14.158(0.079) 

13.337(0.024) 

12.710(0.025) 

8.381(0.024) 

6.086(0.053) 

0.5394E-bO(0.572E-l) 

/ 

146.4630 0.6387 

0.0303 

0.1962 

15.412(0.143) 

14.725(0.200) 

14.232(0.149) 

12.511(0.023) 

12.313(0.024) 

8.115(0.022) 

5.770(0.044) 

0.4890E-(-0(0.127E-bO) 

/ 

146.8068 0.6656 

0.0201 

0.1484 

13.764(0.099) 

12.936(0.125) 

12.150(0.064) 

10.861(0.022) 

10.639(0.020) 

6.588(0.014) 

4.592(0.022) 

0.2862E-i-l(0.192E-b0) 

0.42!lE-i-l(0.174E-H)) 

146.7991 0.7027 

0.0305 

0.1495 

15.785(0.154) 

15.201(0.177) 

14.700(0.167) 

13.204(0.024) 

12.950(0.028) 

9.005(0.033) 

6.099(0.058) 

0.5410E-bO(0.374E-l) 

0.2857E-t-l(0.519E-b0) 

146.7815 0.7380 

0.2618 

0.1485 

16.337(0.172) 

15.473(0.174) 

14.664(0.144) 

14.134(0.042) 

13.850(0.054) 

11.885(0.000) 

8.268(0.000) 


/ 

147.0805 0.7880 

0.2111 

0.1039 

15.820(0.097) 

14.956(0.109) 

14.284(0.097) 

13.734(0.026) 

13.446(0.035) 

11.844(0.000) 

8.758(0.000) 


/ 

149.1699 -0.0233 

0.1392 

0.0299 

15.396(0.090) 

14.640(0.107) 

14.175(0.086) 

13.305(0.025) 

13.108(0.029) 

11.888(0.268) 

8.592(0.000) 


/ 

148.4325-1.0264 

0.1103 

0.0426 

15.413(0.086) 

14.569(0.107) 

13.851(0.074) 

12.767(0.023) 

11.460(0.022) 

8.173(0.022) 

5.372(0.036) 

0.5702E-bO(0.239E-l) 

/ 

148.2377 -0.7920 

0.0898 

0.0334 

14.466(0.073) 

13.847(0.094) 

13.373(0.077) 

12.570(0.023) 

12.507(0.025) 

11.690(0.219) 

8.416(0.000) 


/ 

147.7086 -0.8878 

0.2715 

0.0549 

16.059(0.113) 

15.369(0.145) 

14.596(0.138) 

13.881(0.026) 

13.544(0.032) 

12.440(0.490) 

8.950(0.000) 


/ 

147.4283 -0.8401 

0.0809 

0.0651 

15.132(0.080) 

14.426(0.082) 

13.978(0.076) 

13.298(0.025) 

13.016(0.028) 

11.140(0.131) 

8.647(0.000) 


/ 


Col. 1: coordinates in the epoch of 2000 from Best et al. (2005); Col. 2: redshifts of sources; Col. 3; derived E{B — V)\ Cols. 4—10: magnitudes at J, H, Kg, VVl, W2, W3 and WA bands, respectively, 
with their uncertainty given by the values in parentheses; Cols. 11 and 12; fluxes in units of Jansky at 90 and 140;i.m bands, respectively, with their uncertainty given by the values in pai'entheses. When the value 
in parentheses is zero, the listed magnitudes are the upper limits. We list only the sources in which 90- and 140-/j,m fluxes have high quality such that the source is confirmed and flux is reliable. Only the first 15 
sources are listed here; the full table is available in the electronic version. 


Table 2. Number of detections and upper limits for low-excitation ra¬ 
dio galaxies (LERGs), high-excitation radio galaxies (HERGs) and star¬ 
forming galaxies (SFGs). 


Sub.sample 

Wl(d) 

Wl(u) 

W2(d) 

W2(a) 

M^3(d) 

W3(u) 

V^4(d) 

W4(u) 

LERGs 

1745 

3 

1748 

0 

1032 

716 

220 

1528 

HERGs 

52 

0 

52 

0 

45 

7 

32 

20 

SFGs 

502 

0 

502 

0 

501 

1 

490 

12 


Cols 2—3, 4—5, 6—7 and 8—9 give the number of detections and upper limits for each class at 
four bands of W 1. W2, W3 and W4, respectively. The ‘d’ in parenthesis indicates the number of 
detections, whereas the ‘u’ indicates the number of upper limits. 


the SDSS main galaxy catalogue, which means that objects classi¬ 
fied as quasars by the automated SDSS classification pipeline are 
excluded (Best et al. 2005). 


2.2 Cross-identification 

In recent years, many infrared surveys have been completed by var¬ 
ious instruments at various infrared bands. The Two Micron All 
Sky Survey (2MASS) has swept 99.998% of the whole sky at the 
J (1.25/im), H (1.65/rm) and Ks (2.16/im) bands (Skrutskie et 
al. 2006). The Wide-Field Infrared Survey Explorer (WISE) has 
completed an all-sky survey at the W 1 (3.4/rm), W2 (4.6/rm), W3 
(12/im) and W4 (22/rm) bands and issued two public data releases 
in April 2011 and March 2012. The WISE achieves 5cr point-source 
sensitivities better than 0.08, 0.11, 1 and 6 mjy in the four bands, 
corresponding to Vega magnitudes of 16.5, 15.5, 11.2 and 7.9, re¬ 
spectively; the sensitivity of the WISE W3 (12/rm) band is more 
than 100 times better than that of IRAS (Wright et al. 2010). The 
AKARI (Ishihara et al. 2010), a Japanese infrared space mission, 
has carried MIR to FIR surveys at six bands centred at 9, 18, 65, 
90, 140 and 160 ^m by using two instruments: an infrared cam¬ 
era (Onaka et al. 2007) and a FIR surveyor (FIS; Kawada et al. 
2007). The AKARI achieves higher sensitivity, higher spatial reso¬ 
lution and wider wavelength coverage than the IRAS. 

To obtain the NIR, MIR and FIR photometric data of the stud¬ 
ied radio sources, we primarily used the point source catalogues 
(PSCs) of 2MASS, WISE and AKARI (Skrutskie et al. 2006; Ishi¬ 
hara et al. 2010; Yamamura et al. 2010). The catalogue of the 2712 
RLGs is cross-identified with the 2MASS, WISE and AAAR/-FIS 


PSCs. Most 2MASS PSC sources have spatial resolutions of less 
than ~4 arcsec in all of the three wavelength bands (see 2MASS 
homepag43)- The spatial resolutions of WISE are 6.1, 6.4, 6.5 and 
12.0 arcsec at Wl, W2, W3 and W4, respectively (Wright et 
al. 2010). The spatial resolutions of AKARI are 7 and 48 arcsec 
for the MIR and FIR instruments, respectively (e.g. Matsuta et al. 
2012). We used position offsets of 4, 6 and 20 arcsec to search the 
2MASS, WISE and AKARI-FIS PSC counterparts, respectively, 
for the 2712 RLGs. If more than one object was within the position 
offset radius, the closest object was accepted as the counterpart. 
In practice, only a small number of RLGs fell into the category. 
96% of the 2MASS PSC counterparts and 90% of the WISE PSC 
counterparts were within a 1 arcsec radius of the RLGs; 91% of the 
AAAR/-FIS counterparts were within a 10 arcsec radius. 

Finally, we excluded RLGs without 2MASS counterparts and 
obtained a sample of 2663 RLGs. About 42% of the studied sam¬ 
ples were consistent with a point source. Figure 1 shows the redshift 
distribution of two types of RLGs: radio-loud (RL) AGNs and star¬ 
forming galaxies (SFGs). The classification of RLGs is described 
in section 2.4. The RL AGNs studied are located mainly at redshifts 
of 0.05 < z < 0.3, whereas the SFGs studied are located mainly 
at redshifts of 0.01 < z < 0.15. The infrared photometric data of 
the 2663 RLGs are listed in Table 1, in which Col. 1 gives the coor¬ 
dinates (in the epoch of 2000) from Best et al. (2005); Col. 2 gives 
the redshifts of the sources; Col. 3 gives the derived E(B—V); Cols. 
4—10 give the vega magnitudes at the J,H,Ks,Wl,W2,W3 and 
W4 bands, respectively, with their uncertainties given by the values 
in parentheses (Skrutskie et al. 2006; Ishihara et al. 2010); Cols. 11 
and 12 give the flux values in units of Janskies at the 90 and 140/im 
bands, respectively, with their uncertainties given by the values in 
parentheses (Yamamura et al. 2010). Only the first 15 sources are 
listed here; the full table is available in the electronic version of this 
paper. Table 2 also gives number of detections and upper limits for 
LERGs, HERGs and SFGs. 


2.3 Galactic interstellar extinction correction 

The interstellar extinction corrections caused by our Galaxy, 
E(B — V), are considered in five bands: J, H, Ks, Wl and 

^ http://www.ipac.caltech.edu/2mass/ 











4 


W2. We obtain the values of E{B — V) for all samples through 
the NASA/IPAC extragalactic database (NED Home Pag43), which 
gives the Galactic dust reddening for a line of sight based on 
Schlegel, Finkbeiner and Davis (1998). The extinction of a given 
wavelength band a is usually estimated by A{a) = R{a) x E{B — 
V). Yuan et al. (2013) reported new empirical extinction coeffi¬ 
cients R{a). According to Table 2 of Yuan et al. (2013), we em¬ 
ployed R{a) = 0.72, 0.46, 0.306, 0.19, and 0.15 for the bands 
J, H, Ks, W1 and IV2, respectively, to correct the magnitudes 
of J, B, Ks, Wl and W2. In fact, the extinctions were smaller 
than the uncertainties of the measured magnitudes for the majority 
of the sources. Therefore, we applied interstellar extinction correc¬ 
tions for only about 2 % of our samples. 

2.4 Classiflcation 

At NIR wavelengths, thermal emissions from host galaxies and 
non-thermal emissions from the central nuclei were observed. The 
orientation of AGNs may affect non-thermal emissions from the 
central nuclei. Therefore, to investigate the effect of orientation, 
we used the 13th version of Quasars and AGNs (Veron-Cetty & 
Veron 2010) to search the spectral types of our sample and obtain 
the sources sub-classified as Seyfert 1 (Syl) and Setfert 2 (Sy2) 
galaxies. 

At MIR and FIR wavelengths, thermal emissions from dust 
play a very important role in RLGs. To further investigate the MIR 
properties of different active types of RLGs, we also distinguished 
the active types of RLGs according to Best et al. (2005) and Best 
& Heckman (2012). Based on the location of a galaxy in a plane of 
4000-A break strength (Z)„ (4000)) versus radio luminosity per unit 
stellar mass (Liaghz/MA, the RLGs were divided into RL AGNs 
and SFGs (Best et al. 2005). The radio emissions of RL AGNs are 
mainly generated by jets or lobes, whereas those of SFGs are dom¬ 
inated by star formation. The majority of the samples used by Best 
et al. (2005) are also included in Best & Heckman (2012), although 
a small number of samples were inconsistently classified. For the 
samples overlapping both works, we employed the classification 
from Best & Heckman (2012). Best & Heckman (2012) further di¬ 
vided RL AGNs into HERGs and LERGs; we employed this classi¬ 
fication as well. HERGs and LERGs are believed to belong to dif¬ 
ferent accretion modes, i.e. standard thin discs and ADAFs, respec¬ 
tively (Hardcastle, Evans & Croston 2007). HERGs have higher 
accretion rates than LERGs (Best & Heckman 2012). 


3 NIR PROPERTIES 

Figure 2 shows a (J — H) — {H — Ks) two-colour diagram (top 
panel) and the normalized distributions of the colours H — Ks 
(middle panel) and J — H (bottom panel). In order to easily un¬ 
derstand the spectral shapes of the RLGs, the colours were also 
transformed to their corresponding spectral indices a(Ai, A 2 ). Fol¬ 
lowing Sekiguchi (1987), the spectral index is defined as, 

a(Ai,A2) = log(T’Ai/T’A2)/log(A2/Ai), (1) 

where F\^ and Fx^ are the measured flux densities at wavelengths 
Al and A 2 (A 2 > Ai), respectively. It is noted that the flux densities 
are expressed throughout this paper by using frequency units. The 
corresponding spectral indices of the colours J—H and H — Ks are 

^ http://irsa.ipac.caltech.edu/applications/DUST/ 


also marked in the top panel of Figure 2. In the two-colour diagram, 
a PL ( Ta oc A““) line (dashed line) and a BB line (dotted line) are 
also shown. The PL line can be obtained by calculating the change 
in colour with the spectral index (a). For a PL spectrum, the colour 
in a certain wavelength range AA = A 2 — Ai is derived by 

mAi - mA2 = -2.5 (2) 

where mx^ and mx^ are magnitudes at wavelengths Ai and A 2 , 
respectively, and Toai and ^ 0 X 2 are absolute flux calibrations at 
wavelengths Ai and A 2 , respectively. Similarly, we calculate the 
change in colour with temperature (T) to draw the BB line. For a 
BB spectrum, the colour is derived by 

mAi - mA2 = -2.5 log[(^)(-^)], (3) 

where Bxi and Bx^ are the Planck function with same temperature 
at wavelengths Ai and A 2 , respectively. 

In addition, following Kouzuma & Yamaoka (2010), the top 
panel of Figure 2 shows the intrinsic loci of dwarf stars (magenta 
line), giant stars (red line) and Classical T Tauri Stars (CTTS; blue 
line) according to Bessell & Brett (1988) and Meyer et al. (1997). 
The reddening vector from Rieke & Lebofsky (1985), which in¬ 
dicates the shifting direction of the stellar and CTTS loci on the 
(J — H) — {H — Ks) diagram, is also plotted. The infrared colour 
of the RLGs is different from the colours of various types of stars. 
Kouzuma & Yamaoka (2010) proposed the following NIR colour 
criteria for selecting AGNs: 

(J - 77) < 1.7(77-iT,)-0.118, (4) 

(J - 77) < 0.61(77 - A",)-P 0.5. (5) 

Equations 4 and 5 are also indicated in the two-colour diagram by 
dot—dashed and double-dot-dashed lines, respectively. 64% of RL 
AGNs and 48% of SFGs achieved these selection criteria, which 
is a lower percentage than that given by Kouzuma & Yamaoka 
(2010) for low-redshift AGNs. It is noted that the sample used in 
the present study was located at redshifts of 0.01 < z <0.56 where 
the NIR radiation from the centre of the AGN should be contam¬ 
inated by host galaxies with intrinsic colours similar to those of 
normal stars. On the contrary, the majority of the samples stud¬ 
ied by Kouzuma & Yamaoka were quasars that had more luminous 
cores than those of our sample. Hence, the NIR selection was less 
effective for the RL AGNs than for other types of AGNs. 

It should be noted that most of the sources achieving the selec¬ 
tion criteria were distributed in the region around the PL line, which 
indicates that their NIR spectra can be approximately described by 
a PL spectrum. The top panel of Figure 2 also shows that the major¬ 
ity of the RL AGNs and SFGs had indices of —2 < q:(77, Ks) < 1 
and —2 < a( J, 77) < 0.4, which is less than the value of 4 sup¬ 
posed by the Rayleigh—Jeans spectra. AGNs produced emissions 
at nearly all wavelengths, and their emissions were dominated by 
non-thermal radiation. The ultraviolet to MIR continuum of AGNs 
is represented by a PL (Stem et al. 2005). However, host galaxies 
may contaminate the NIR colour of AGNs. The spectra from the 
star population in the host galaxy were determined by using the 
BB spectra and peak at approximately 1.6 /rm (Elvis et al. 1994; 
Stem et al. 2005; Richards et al. 2006; Park et al. 2010). For exam¬ 
ple, Dultzin-Hacyan & Benitz (1994) pointed out that the dominant 
NIR emissions of Sy2 galaxies are primarily produced by evolved 
massive red giants and supergiants in the host galaxies; De Breuck 
et al. (2010) used the spectrum of a stellar population as a com¬ 
ponent in fitting the NIR and MIR spectra of high-redshift radio 
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Figure 2. (J — H)—(H — it's) two-colour diagram (top panel) and the nor¬ 
malized distribution of the colour H — Ks (middle panel) and J-Ft (bottom 
panel). The power-law spectral indices a{J, H) and q(H,Ks) are also 
marked in the top panel. A BB line with temperatures and a PL line with 
spectral indices are plotted as the dotted and dashed line, respectively. In 
the middle and bottom panels, black and red lines represent RL AGNs and 
SFGs, respectively. 


galaxies. To study the degree in which objects with non-thermal 
emissions and stellar populations affect the NIR colour of active 
galaxies, Kouzuma & Yamaoka (2010) first examined the distribu¬ 
tion of four types of objects in the (J — H) — {H — Ks) dia¬ 
gram including microquasar candidates, low-mass X-ray binaries, 
cataclysmic variables and massive young stellar objects and found 
that most of the objects were not close to the PL line and did not 
achieve the NIR colour criteria of AGNs. Based on these results, 
they determined that the dominant NIR radiation of the four ob¬ 
jects was likely thermal emission. They further examined the loci 
of normal galaxies without considering the contribution from the 
nuclei in the (J — H) — {H — Ks) diagram and concluded that the 
normal galaxies with 0.6 < z < 0.8 have intrinsic colours similar 
to those of normal stars, whereas galaxies with z ~ 1 achieve the 
NIR colour criteria of AGNs. Therefore, the probable objects with 
non-thermal radiation and the stellar populations in the host galax¬ 
ies cannot cause active galaxies with z < 0.8 to approach the PL 
line. Thus, if active galaxies with z < 0.8 are close to the PL line 
in the (J — H) — [H — Ks) diagram, their NIR emissions should 
be dominated by non-thermal emissions associated with the nuclei. 
For the SFGs achieving the selection criteria (48% of SFGs), the 
non-thermal emissions from the nuclei dominate their NIR emis¬ 
sions. It should also be noted that a great number of sources are 
distributed in the region around the BB line, which indicates that 
their NIR emissions may be dominated by thermal radiation be¬ 
cause the BB line is considered as an indicator of thermal radiation 
from the starlight in the host galaxy. The contribution ratio of the 
centre nuclei and their host galaxies to NIR emissions determines 
the position of active galaxies in the (J — H) — {H — Ks) diagram. 

Chen & Shan (2009) also reported that almost all type 1 
quasars are distributed in the region around the PL line in the 
(J — H) — (H — Ks) diagram due to their predominantly non- 
thermal properties, whereas the majority of Sy2 galaxies are dis¬ 
tributed in the region around the BB line due to the obscuring of 
the non-thermal component of AGNs by the dust torus. Here the 
(J — H) — (H — Ks) diagram shows that many Syl RLGs tend 
to lie close to the BB line. The difference in distribution of Syl 
RLGs and type 1 quasars occurs because that type 1 quasars have 
more luminous nuclei; therefore, the NIR emissions from the nuclei 
are significantly stronger than those from their host galaxies. It is 
noted that a great number of Sy2 RLGs are distributed in the region 
around the PL line. A similar case also exists for type 2 quasars 
(Chen & Shan 2009). Chen & Shan (2009) determined that about 
43% of type 2 quasars are distributed in the region around the BB 
line, whereas about 57% of type 2 quasars are distributed near the 
PL line. They cited the possibility that some NIR emissions from 
the central nuclei can penetrate the dust torus because it is not sub¬ 
stantially thick and the viewing angle is not very large. Another 
possible reason is that the torus is clumpy; therefore, the radiation 
from the central nuclei can penetrate the clumpy torus at a certain 
angle. 

The middle and bottom panels of Figure 2 show the normal¬ 
ized histograms of the colours J — H and H — Ks of RL AGNs 
and SFGs, respectively. In addition, we performed a Kolmogorov- 
Smirnov (K-S) test to evaluate the similarity in the distributions of 
these colours for RL AGNs and SFGs. We determined that the sim¬ 
ilarity of J — H has a significance level of less than 5%, whereas 
that of H — Ks is larger than 5%. These results indicate that the 
distributions of J — 77 for RL AGNs and SFGs are essentially dif¬ 
ferent, whereas those of H — Ks are essentially the same. Thus, a 
systematic difference exists in J-band emissions of RL AGNs and 
SFGs, the reason for which needs to be investigated. 
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Figure 3. ([3.4]-[4.6])-([4.6]-[12]) (top panel) and ([12]-[22])- 
([4.6]—[12]) (bottom panel) two-colour diagrams. A PL line marked with 
spectral indices and a BB line are also plotted by dashed and dotted lines, 
respectively. The upper limits are indicated as arrows. 


4 MIR COLOUR PROPERTIES 
4.1 Colour properties 

Figure 3 shows ([3.4] - [4.6]) - ([4.6] - [12]) and ([12] - [22]) - 
([4.6] — [12]) two-colour diagrams. In each diagram, a PL line and 
a BB line are plotted with a dashed and dotted line, respectively. 

From Figure 3, we can see that RL AGNs and SFGs are clus¬ 
tered in distinctly separate regions. Compared to the RL AGNs, 
the SFGs have a redder [4.6] — [12] colour and a slightly red¬ 
der [3.4]-[4.6] colour. The PL spectral index a(3.4, 4.6) of the 
RL AGNs changes from approximately 0 to 2, q( 4.6, 12) changes 
from approximately —2 to 2 and a?( 12 , 22 ) changes from approx¬ 
imately —4 to 0. This suggests that their observed spectra behave 
differently within the wavelength range of 4.6—22/rm. One type is 
q( 4.6, 12) < 0 and q( 12,22) < 0, showing an increase of F\ 
with A from 4.6 to 22 /rm, and the other type is q:( 4.6, 12) > 0 and 
a(12,22) < 0, showing a ‘turning point’ in the spectra, i.e. F\ 
decreases with A from 3.4 to while F\ increases with A from 


a(4.6,12) 



Figure 4. Probability distribution of LERGs and SFGs. Grey-scale and 
30 isodensity contour lines (black solid lines: LERGs; black dotted lines: 
SFGs) indicate the probability distribution of LERGs and SFGs. HERGs 
are superimposed as red circles in this figure for comparison. A PL line 
marked with different spectral indices and a BB line are also plotted with 
dashed and dotted lines, respectively. Blue solid lines indicate the loci of 
models with T = SbOOA and a = —2, —3 and —4. Red solid lines 
indicate the loci of models with T = bSOOFf and a = —2, —3 and —4. 
Numbers along the loci are the fraction of contribution of the PL component 
to the VF2-band flux. Objects with upper limits are excluded. 

12 to 22/rm. However, the SFGs show red spectra, i.e. Fx increases 
with A from 4.6 to 22/rm. 

To understand the MIR colours of the RLGs, a BB component 
and a PL component are used. The BB component roughly models 
the component of the stellar population, whereas the PL compo¬ 
nent represents the contributions of multi-temperature hot dusts or 
non-thermal emissions associated with the centre nuclei. The ultra¬ 
violet to MIR continuum of AGNs is described by a PL (Stern et 
al. 2005). The thermal emissions of multi-temperature hot dusts are 
also described by a PL, as seen in Figure 3 of Brandi et al. (2006). 
For example. Smith et al. (2007) used the modified blackbodies at 
fixed temperatures (35, 40, 50, 65, 90, 135, 200 and 300 K) as ther¬ 
mal dust continuum components to decompose the MIR spectra of 
SFGs. The two components are composed according to the fraction 
of their contribution to the VF 2 -band flux, so we can calculate the 
flux densities of the composed model at different wavelengths by 

Fx = Bx{T) + AixX-°‘. (6) 

Here, Bx(T) represents the Planck function of a BB with temper¬ 
ature T and Ai = -r^i-where 71 is the ratio of the 

1 —n A7? „ ' 

radiation flux density of the PL component to the total flux den¬ 
sity at the W 2 band. For given values of temperature T and spec¬ 
tral index a, we can obtain the dependence of the colour of the 
composed model on 71 and plot the loci of the composed model 
with blue and red lines, as in Figure 4. Here, for simplicity, we 
assign the BB temperatures T = 3500 and 5500 K, according to 
the later stellar effective temperature. We use different spectral in¬ 
dices a to model the potential energy distributions of emissions 
from multi-temperature dusts or non-thermal emissions associated 
with the centre nuclei. From the ([12]-[22])—([4.6]-[12]) diagram, 
we can see that the spectral indices a( 12 , 22 ) of few sources are 











Fraction Fraction Fraction 


Infrared colour properties of nearby radio-luminous galaxies 1 




[4.6H12] 


0.25 


--- 


. LERGs1 





- HERGsl 

0.20 

1 _' 1 



0.15 

-I'" i 



0.10 




0.05 




_ 1 _ , _ 1 _J _ i_ L, _ 1 _,_ 1 _ 

1 1 ' 

1 1 1 


1.0 1.5 2.0 2.5 3.0 3.5 4.0 

[12]-[22] 

Figure 5. Normalized histograms of three colours of HERGs (red) and 
LERGs (black). Top panel is for colour [3.4] — [4.5], middle panel is for 
colour [4.6] — [12], and bottom panel is for colour [12] — [22]. The objects 
with upper limits are excluded. 


less than -4, so we set a ^ — 4 , such as a = —2, —3, and — 4 , 
respectively. Numbers along the loci are the fraction of the contri¬ 
bution of the PL component to the VK2-band flux. Figure 4 also 
shows the probability distribution of LERGs and SFGs with grey¬ 
scale and 30 isodensity contour lines (black solid lines for LERGs 
and black dotted lines for SFGs). The probability distribution is 
evaluated using the kernel density estimation technique, which is 
an effective method of estimating the probability function of a mul¬ 
tivariate variable without any assumption. For comparison, we also 
superimpose HERGs as red circles in Eigure 4. From Figure 4, we 
can see that: the smaller the spectral index a, the redder the colour 
[4.6]-[12]; the larger the contribution of the PL component to the 
W2-band flux, the closer to the PL line the distribution of the RLGs 
in the ([4.6]-[12])-([12]-[22]) diagram is. Although our parameters 
cannot cover the distribution of all the RLGs, the composed model 
helps us to understand the role of the BB and PL components in the 
MIR colours of RLGs. 

Figure 4 also shows that LERGs have a distribution of a 
double-core structure, which means that LERGs do not have singu¬ 
lar infrared properties. The main difference between the two cores 
is their different [ 3 . 4 ] — [ 4 . 6 ] colour. More spectral data are required 
to understand this difference. It is noted that the W 1 band covers 
the PAH emission at 3.3 fim and the W3 band covers the 9.8-/rm 
silicate absorption and the 1 L3-/rm PAH emission, while there are 
no PAH emissions in the W2 and W4 bands (e.g. Brandi et al. 
2006 ; Smith et al. 2007 ). Therefore, many LERGs display a bluer 
[ 3 . 4 ] — [ 4 . 6 ] colour than the red line in Figure 4, due to the 3.3-/rm 
PAH emission. 

Similarly composed models are used to mimic the distribu¬ 
tion of most RL AGNs and SFGs in the ([12] - [22]) - ([22] - 
[12]) diagram. Although the ([ 3 . 4 ] — [ 4 . 6 ]) — ([ 4 . 6 ] — [12]) and 
([12] — [22]) — ([22] — [12]) diagrams cannot definitively confirm 
the nature of thermal or non-thermal emissions of RL AGNs and 
SFGs at MIR wavelengths, they help us understand the spectral 
shapes of RL AGNs and SFGs. For SFGs, we presume that multi¬ 
temperature hot dusts dominate MIR emissions, as what Brandi et 
al. (2006) pointed out. The ([12] - [22]) - ([22] - [12]) diagram 
also shows that the [12] — [22] colour of RL AGNs is more widely 
distributed than that of the SFGs, which suggests that the MIR radi¬ 
ation components of RL AGNs is more diverse than that of SFGs. 
MIR radiation of RL AGNs may include the non-thermal emis¬ 
sions associated with AGNs and the thermal emissions from torus 
or host galaxies. The two responsible factors need be further quan¬ 
titatively studied for RL AGNs. Stern et al. (2012) and Assef et al. 
(2013) presented a simple MIR colour criterion of selecting AGNs 
[ 3 . 4 ] — [ 4 . 6 ] > 0.8. In our studied samples, 0.06% of RL ANGs (14 
RL AGNs) and 10% of SFGs (50 SFGs) achieve this criterion. 

In addition, from the ([ 3 . 4 ] — [ 4 . 6 ]) — ([ 4 . 6 ] — [12]) and 
([12] — [22]) — ([ 4 . 6 ] — [12]) diagrams, we can see the obvious dif¬ 
ference in the distribution of LERGs and HERGs, though LERGs 
and HERGs partly overlap. To show the difference, we plot normal¬ 
ized histograms of three MIR colours ([ 3 . 4 ] — [ 4 . 6 ], [ 4 . 6 ] — [12], 
and [12] — [22]) of LERGs and HERGs in Figure 5. This plot shows 
that HERGs have redder [ 3 . 4 ] — [ 4 . 6 ] and [ 4 . 6 ] — [12] colours on 
average. LERGs are mainly located in the region of [ 3 . 4 ] — [ 4 . 6 ] 
< 0.4 and [ 4 . 6 ] — [12] < 0.8, whereas a large number of HERGs 
have [ 3 . 4 ] — [ 4 . 6 ] > 0.4 and [ 4 . 6 ] — [12] > 0.8. In particular, few 
HERGs have [ 3 . 4 ] - [ 4 . 6 ] < 0.1 and [ 4 . 6 ] - [12] < 0.8. HERGs 
have a narrower colour [12]-[22] distribution than LERGs. Eigures 
3 and 4 show that HERGs are distributed in the region where MIR 
emissions are dominated by the power-law component. Eigures 4 
also shows that a large number of LERGs are distributed in the re- 
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Figure 6. Plots of the logarithms of 4.6-fim luminosity versus the colour 
[3.4] — [4.6] (top panel), the logarithms of 12-/rm luminosity versus the 
colour [4.6] — [12] (middle panel), and the logarithms of 22-fim luminosity 
versus the colour [12] — [22] (middle panel). 


gion dominated by the steller population and a small number of 
LERGs have redder [4.6] — [12] colour. The difference between 
LERGs and HERGs may be concerned with their different accre¬ 
tion modes. It is suggested that HERGs and LERGs have differ¬ 
ent accretion modes (Hardcastle et al. 2007). LERGs are suggested 
to be fuelled via ADAEs (hot-mode accretion) with low accretion 
rates, whereas HERGs are associated with a standard thin accretion 
disc (cold-mode accretion) at high accretion rates (Shakura & Sun- 


Table 3. Results of partial correlation analyses between luminosities and 
colours. 


Abscissa 

Ordinate 

Subsample 

Number 

P.z 


[3.4] - [4.6] 


LERGs 

1748 

0.459 

8.80 



HERGs 

52 

0.585 

6.50 



SFGs 

502 

0.460 

11.94 

[4.6] - [12] 


LERGs 

1748 

0.284 

27.80 



HERGs 

52 

0.617 

13.0 



SFGs 

502 

0.197 

11.39 

[12] - [22] 

logic (^^22^xm) 

LERGs 

1748 

0.027 

7.90 



HERGs 

52 

0.179 

3.35 



SFGs 

502 

0.388 

17.29 


Notes. p,z is Kendall’s partial correlation coefficient with its associated null probability, t j g '> Z 
means that the correlation is significant at the 3cr level. 


yaev 1973; Hardcastle et al. 2007; Best & Heckman 2012). HERGs 
with high accretion rate have luminous core, while LERGs with 
low accretion rate have dim core. Therefore, the MIR emissions of 
HERGs are dominated by non-thermal emissions associated with 
AGNs; the MIR emissions from central core of LERGs are signifi¬ 
cantly contaminated by stars and dusts of host galaxies. Accretion 
rates play an important role in the switch between hot- and cold¬ 
mode accretions (Esin et al. 1997). However, accretion rates are 
also modulated by feedback from central AGNs, and different ac¬ 
cretion modes of the central AGNs have different effects on this 
feedback (Ciotti & Ostriker 2001; Gan et al. 2014). Determining 
what the key of switch is between hot- and cold-mode accretions is 
still an open issue. 


4.2 Luminosity versus colour 

Here, we check possible correlations between MIR luminosities 
and colours for different types of RLGs. We calculate the lumi¬ 
nosity distance (79l) by using the redshift listed in Table 1 (cos¬ 
mology parameters are assumed in Section 1) and then obtain the 
absolute luminosities. The K correction may become important for 
high-redshift objects. We apply the K correction to our sample. Eor 
simplicity, MIR spectra of our sample are approximated by sim¬ 
ple PL form (oc The intrinsic luminosities L{u) are calculated 
as follows: L{u) = 4:TvDlf{u){l -|- (e.g. Bourne et al. 

2011 ), where /(t^) is the observed flux density at the observing fre¬ 
quency iz. Here, a can be calculated by using Equation 1. In Eigure 
6 , we plot the logarithms of 4.6-fim, and 22-fim luminosi¬ 

ties versus the [3.4] — [4.6], ]4.6] — [12], and [12] — [22] colours, 
respectively. We can see that the 12-/rm and 22-fim luminosities 
of HERGs are on average higher than those of LERGs. Similar re¬ 
sults are also shown in Eigure 7 of Giirkan et al. (2014). Giirkan et 
al. (2014) remarked that a division of LERGs and HERGs clearly 
stands out in the 22-^m luminosity. However, based on our studied 
sample, HERGs cannot be separated from LERGs using the 22-fim 
luminosity. 

In order to test physical correlations between luminosities and 
colours, we further calculate Kendall’s partial correlation coeffi¬ 
cient (p.z) with censored data using the methods by Akritas & 
Siebert (1996). p.z measures the degree of association between 
two random variables, with the effect of a set of controlling ran¬ 
dom variables removed, p.z is especially important for excluding 
the effect of redshift. The values of p.z are within —1 and 1. The 
larger their absolute values, the stronger the correlations. The pos¬ 
itive values indicate the positive correlation and the negative val¬ 
ues indicate the negative correlation. Table 3 gives results of partial 
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Figure 7. Histogram of colour [4.6] — [12] of RL AGNs (black dotted line) 
and SFGs (red dashed line). The objects with upper limits are excluded. 


correlation analyses. For LERGs, a physical relationship presents 
between logj^Q(AL4,6/im) and [3.4] — [4.6] and the correlation of lu¬ 
minosities vs. colours becomes weak with increasing wavelength. 
As shown in the bottom panel of Figure 6, no correlation between 
logj^Q(AZ/22Mm) versus [12] — [22] is seen for LERGs. For HERGs, 
good physical correlations exist between log^Q(AI/4.6^m) versus 
[3.4] — [4.6] and between logj^Q(ALi2^m) versus [4.6] — [12], but 
the correlation of logj^Q(AL22(jm) versus [12] — [22] is weaker. For 
SFGs, weak correlations present between logj^g(AL4,6/im) versus 
[3.4] — [4.6] and between logj^g(AL22^m) versus [12] — [22]. 


4.3 Distinguishing RL AGNs and SFGs by MIR colonrs 

The activity in galaxies is thought to be influenced by two compo¬ 
nents, star formation and AGNs. Some galaxies are dominated by 
either star formation or the AGN, while others exhibit comparable 
contributions from both components. The AGN of NGC 1068 is 
surrounded by a ring of star formation that contributes half of the 
total luminosity of NGC 1068 (Telesco & Decher, 1988). The two 
components are important in the understanding of galaxy energet¬ 
ics and galaxy evolution. It is necessary to distinguish the dominant 
component of the activity for active galaxies. Both starburst and 
Sy2 AGN galaxies exhibit emission lines with FWHM< 1000 km 
s“^. However, Baldwin et al. (1981), Veilleux & Osterbrock (1987) 
and Kauffmann et al. (2003) demonstrated that Sy2 AGN galaxies 
may be distinguished from normal SFGs based on the diagnostic 
diagram of [O III]/H/3 versus [HII]/Ha, though not all galaxies can 
be easily segregated. Machalski & Condon (1999) separated the RL 
AGN and star-forming populations on the basis of FIR-radio flux 
ratios, FIR spectral indices and radio-optical flux ratios. Best et al. 
(2005) pointed out that RL AGNs are fairly well separated from 
SFGs in the plane of Z)„(4000) versus Liaghz/M t. 

Based on the classification provided by Best et al. (2005) and 
Best & Heckman (2012), we plot a histogram of the [4.6] — [12] 
colour in Figure 7. Figure 7 shows the separation of RL AGNs from 
SFGs; 95% of RL AGNs have [4.6] - [12] < 3.0, and 94% of SFGs 
have [4.6] - [12] > 3.0. The colours [3.4] - [12], [3.4] - [22], 
and [4.6] — [22] have similar functions in the separation of the two 
types of galaxies, which is attributed to the fact that the 12-/rm and 
22-^m emissions of SFGs are redder than those of RL AGNs. 


5 FIR COLOUR PROPERTIES 

We search AKARI-FIS PSC and obtain AKARI-FIS data with high 
quality. AAAR/-FIS data have four-level flux quality indicator. High 
quality, i.e. the flux quality indicator equals 3, means that the source 
is confirmed and the flux is reliable. We neglect the data with the 
flux quality indicator equal to 1 and 2, which means that the source 
is not confirmed or the flux is not reliable. 15 and 282 AKARI 90- 
^m counterparts are obtained for RL AGNs and SFGs, respectively. 
The detection rate of RL AGNs at the 90-^m band is ~ 0.7%, 
which is much less than ~ 56% of SFGs. The majority of the 
RL AGNs studied here are not included in the AAARZ-FIS PSC, 
though the AAAR/-FIS survey covered more than 94% of the entire 
sky, and the AKARI-FIS PSC includes 0.4 million infrared sources 
brighter than the detection limit (approximately 0.6 Jy) at 90 ^m 
(Yamamura et al. 2010). This suggests that the 90 ^m flux density 
of most of RL AGNs is lower than 0.6 Jy. In addition, 15 AKARI 90- 
/rm counterparts of RL AGNs include nine LERGs and one HERG, 
respectively. The detection rate of LERGs and HERGs reach ~0.5 
and ~2 per cent at 90-/im band, respectively. The FIR emissions 
at 90 and 140 /rm are thought be associated with cold dust in host 
galaxies (e.g. Netzer et al. 2007; Mullaney et al. 2011), and then 
HERGs tend to have more abundant cold dust than LERGs, on av¬ 
erage. 

The a(22,90) - a(12,22) and a(90,140) - a(22, 90) di¬ 
agrams are also shown in Figure 8 with PL and BB lines plotted 
as dashed and dotted lines, respectively. A set of solid lines in the 
a(90,140) — a(22, 90) diagram indicate the loci of the composed 
model, with a cold disc component of ~30 K and a warm dust 
component of 80—100 K. The two components are composed ac¬ 
cording to the fractional contribution of the warm dust component 
to the 90 ^m-band flux. We can calculated the flux densities at a 
given wavelength by 

Fx = Bx{Tc) -b Az X Bx{T^). (7) 

Here, B\{Tc) and B\{T,i) represent the Planck functions of the 
cold component with temperature B and the warm component with 
Tw, respectively, and Az = , , where 7z is the ratio 

of the radiation flux density of the warm component to the total 
flux density at 90 fim. We can then draw the set of black line in the 
a(90,140) — a(22, 90) diagram by computing the spectral index 
between two bands from Equation 1 when 7z changes. 

Figure 8 shows a'(22, 90) < 0 for all RL AGNs detected at 
90jrm, which indicates that some RL AGNs have a red FIR colour. 
In combination with Figure 3, we can see that the spectral in¬ 
dices for the majority of SFGs are q:(3.4, 4.6) >0, q( 4.6, 12) <0, 
a(12, 22) <0, a(22, 90) <0 , and a(90,140) <0, i.e. their spec¬ 
tral indices do not cross zero, which suggests that SFGs have sim¬ 
ilar spectra from 3.4 to 90/im. The spectra may be described by a 
PL with a similar spectral index, as shown by the 5—38/rm MIR 
spectra of 22 starburst nuclei presented by Brandi et al. (2006, see 
Figure 3 in their paper). However, the bottom panel in Figure 8, 
the a(90,140) — a(22, 90) diagram, shows that a small number 
of SFGs with a(90,140) < —1.4 are distributed below the black 
solid lines and close to the PL line. A set of black solid lines, 
plotted in the a(90,140) — a(22, 90) diagram, indicates the com¬ 
posed model of a cold disc component (~30 K) and a warm dust 
component. As shown in the a(90,140) — a(22, 90) diagram, the 
sources above the black solid lines and close to the BB line may 
be explained by thermal emissions from dust at MIR wavelengths, 
whereas the MIR properties of the source distributed in the region 
around the PL line cannot be explained by thermal emissions of 
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Figure 8. The o(22, 90) - o(12, 22) and o(90,140) - o(22, 90) dia¬ 
grams. A BB line marked with temperature and a PL line are also plot¬ 
ted as the dashed and dotted line, respectively. A set of solid lines in the 
0(90,140) — a(22, 90) diagrams indicate the loci of models with a cold 
dust component of 30 K and a warm dust component of 80—100 K. Circles 
represent RL AGNs and triangles represent SFGs. 

the composed model, meaning that non-thermal emissions may be 
used to understand their MIR properties. In general, the FIR emis¬ 
sions of SFGs are thought to be caused by the heating of dust by 
young, massive stars (Condon et al. 1991; Helou & Bicay 1993; 
Machalski & Condon 1999). The q:(90, 140) — q;( 22, 90) diagram 
suggests that there are many SFGs whose FIR emissions may not 
come from dust. We also suggest the criteria q( 90, 140) < —1.4 
to select active galaxies with non-thermal emissions at MIR wave¬ 
lengths. 


6 SUMMARY 

The radio emissions of RLGs are dominated either by AGNs or 
by supernova remnants in their star-forming regions; however, both 
the dusty tori associated with the AGNs and the dust and gas clouds 
heated by young stars emit MIR or FIR emissions. Here, we analyse 
the infrared colour properties of a large sample of RLGs based on 
the 2MASS, WISE and AKARI survey data. RL AGNs and SFGs 
are separately distributed in a ([3.4] — [4.6]) — ([4.6] — [12]) dia¬ 
gram, which provides the colour [4.6] — [12] to separate RL AGNs 
from SFGs. 95% of RL AGNs have [4.6] - [12] < 3.0, and 94% 


of SFGs have ]4.6] — [12] >3.0. This difference may be due to the 
signihcant contribution of the dusts of SFGs to MIR emission. It is 
noted that a small number of RL AGNs are distributed in the region 
of SFGs in the ([3.4] — [4.6]) — ([4.6] — [12]) diagram, which in¬ 
dicates that the RL AGNs in the region of SFGs have redder MIR 
and FIR colours than ordinary RL AGNs. The 3.4-90^m observed 
spectra of SFGs are described by similar spectra, like a PL, which 
are indicated by the distribution of SFGs around the PL line in 
(]12] - ]22]) - ([4.6] - [12]) and a(90,140) -a(22, 90). However, 
the 3.4—22/rm spectra of RL AGNs shows different types, due to 
the variation ofQ:(4.6,12) from —2 to 2. One type has the flux den¬ 
sity increasing with wavelength from 3.4 to 22/im; the other type 
has the flux density decreasing with wavelength from 3.4 to 12/rm, 
whereas flux density increases with wavelength from 12 to 22/im. 

The MIR colours of RLGs are described by a model with a 
BB component and a PL component. However, the PL component 
cannot exclude the contribution from thermal emissions, because 
the PL component may be described by multi-colour BB spectra, 
indicating dusts with various temperatures. 

LERGs and HERGs have different distributions of the [3.4] — 
]4.6], [4.6] — ]12] and [12] — [22] colours, which means that LERGs 
and HERGs have intrinsically different MIR properties. 

The a(90,140) — q( 22,90) diagram indicates that many 
SFGs with a (90,140) < —1.4 are distributed in the region around 
the PL line. They cannot be described by the composed model 
with a cold disc component (~30 K) and a warm dust compo¬ 
nent involved in star formation. The FIR emissions of the SFGs 
with q:( 90, 140) < —1.4 may not originate from dusts but instead 
from non-thermal processes. Therefore, the criterion a (90,140) < 
— 1.4 provides a way to select active galaxies, for which non- 
thermal emissions may contribute to FIR wavelength. 
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